Track estimation with binary derivative 

observations 

Adrien Ickowicz 
Abstract 

We focus in this paper in the estimation of a target trajectory defined by whether a time constant parameter 
in a simple stochastic process or a random walk with binary observations. The binary observation comes from 
binary derivative sensors, that is, the target is getting closer or moving away. Such a binary obervation has 
a time property that will be used to ensure the quality of a max-likelihood estimation, through single index 
model or classification for the constant velocity movement. In the second part of this paper we present a 
new algorithm for target tracking within a binary sensor network when the target trajectory is assumed to be 
modeled by a random walk. For a given target, this algorithm provides an estimation of its velocity and its 
position. The greatest improvements are made through a position correction and velocity analysis. 

I. Introduction 

Sensor networks are systems made of many small and simple sensors deployed over an area in an attempt 
to sense events of interest within that particular area. In general, the sensors have limited capacities in 
terms of say range, precision, etc. The ultimate information level for a sensor is a binary one, refeiTing to 
its output. However, it is important to make a distinction according to the nature of this binary information. 
Actually, it can be related to a — 1 information (non-detection or detection) or to relative {— , +} motion 
information. For example, if the sensors are getting sound levels, instead of using the real sound level 
(which may cause confusion between loud near objects and quieter close objects), the sensor may simply 
report whether the Doppler frequency is suddenly changing, which can be easily translated in whether 
the target is getting closer or moving away. Moreover, low-power sensors with limited computation and 
communication capabilities can only perform binary detection. We could also cite video sensors, with the 
intuitive reasoning: the target is getting closer if its size is increasing. The need to use that kind of sensor 
networks leads to the development of a model for target tracking in binary sensor networks. 

We consider a sensor network, made with N sensors (e.g. video),with (known) positions. Each sensor can 
only gives us a binary { — , +} information f2l, i.e. whether the target-sensor distance is decreasing (— ) or 
increasing (+). This "choice" can result from severe communication requirements or from the difficulties 
from fusing inhomogeneous data. Even if many important works deal with proximity sensors Q, (|6l, we 
decide here to focus on the binary {— ,+} information |2l. Here, the aim is to estimate the parameters 
defining the target trajectory. Even if our methods can be rather easily extended to more complex models 
of target motion, we decide to focus here on a constant velocity movement. Actually, this framework is 
sufficiently general to present the main problems we have to face, as well as the foundations of the methods 
we have to develop for dealing with these binary data. See fig. Ufor an example. 

In a first time, the observability requirements are considered. Then, we turn toward the development of 
specific estimation methods. Especially, the new concept of the velocity plane is introduced as an exhaustive 
representation of the spatio-temporal sequence of binary data. It is then used both in a separation-oriented 
framework (SVM) and in a projection pursuit regression (PPR) one. The corresponding methods are carefully 
presented and analyzed. 

In the following pait we release the assumption of (piecewise) constant velocity motion, and we try to 
follow both position and velocity in real time. In paiticular, it is shown that it is the trajectory "diversity" 
which renders this possible. 

Obviously, tracking a diffusive Markovian target widely differs from the (batch) estimation of deterministic 
parameters. However, both problems present strong similarities. Indeed, the geometrical propeities remains 
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Fig. 1 . A view of a getting closer car 



the same at each instant. Once the target motion model has been introduced, the most important properties we 
used to perform the tracking are presented. Then, the method which allows us to perform adapted corrections 
for tracking the target is presented. It is the main contribution of this part of the paper. 

Simulation results illustrate the behavior of the estimators, as well as the performances of the tracking 
algorithm. We conclude on further works about the tracking in binary sensor networks. 



II. Binary Sensor Network Observability Properties 

Let us denote S; a sensor whose position is represented by the the vector Similarly, the vector xt 
represents the position vector of the target at the time-period t. Let us denote di{t) the (time-varying) 
distance from sensor S; to the target at time t. Then, we have that: 

d,{t) \^ d^{t) <0 ,or. (xt-t,,vt) <0, (1) 

where Vf is the instantaneous target velocity. We thus have the following lemma. 

Lemma 1. Let Si (resp. Sj) a sensor whose the target distance is decreasing (resp. increasing) at the 
time-period t, then we have: 

(tj,vt) < (x,,vt) < (t„vt) . (2) 



If we restrict to binary motion information, we consider that the output Si{t) of a sensor (at time t) is +1 
or —1 according to the distance di{t) is decreasing or increasing, so that we have: 

s,(t) = +1 if d^{t) < , 

(3) 

Sj{t) = -1 if dj{t) > . 

Let us denote A the subset of sensor whose output is +1 and B the subset of sensors whose output is — 1, 
i.e. A = {si\ Si{t) = +1} and B = {sj\ Sj{t) = —1} and C{A) and C{B) their convex hulls, then we have 
0: 

Proposition 1. C{A) H C{B) = and xj <^ C{A) [J C{B). 

Proof: The proof is quite simple is reproduced here only for the sake of completeness. First assume that 
C{A) Pi C{B) 7^ 0, this means that there exists an element of C{B), lying in C{A). Let s be this element 
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(and t its associated position), then we have (t <£ C{B)): 



t = t, , > and J2(3, - 1 

o that we have on the first hand: 

t,vt) = (tj,vt) < (xf,vt) (seeeq.EI, 

jeB 

ind, on the other one {t e C{A)y. 

= (tj,Vf) > ( j 



j&B jGB 

SO that we have on the first hand: 

rj.] ( see en. FtI). 
and, on the other one {t e C{A)y. 

{(t„v(t))} > (Xt,Vt) 



Thus a contradiction which shows that C{A) f] C{B) = 0. For the second part, we have simply to assume 
that x(t) e C(A) ( Xf = ti, ttj > 0), which yields: 

(xt,vt) = y^Oi (ti,Vf) > min(tj,vt), (5) 

— ^ zGA 

which is clearly a contradiction, idem if X{t) E C{B). 
□□□ 

So, C{A) and C(i3) being two disjoint convex subsets, we know that there exists an hyperplane (here a 
line) separating them. Then, let Sfc be a generic sensor, we can write tfc = A + /i v^, so that: 

(tfc,vt) =A||vt||' >0^ A>0. (6) 

This means that the line spanned by the vector separates C{A) and C{B). Without considering the 
translation and considering again the {vf,v^} basis , we have : 

tfe e A llvtll^^ > (xt,vt) , 



tfc e A llvtir < (xt,vt) . 

Thus in the basis (vf,v^), the line passing by the point | ^-^ — s^jO ) ™d whose direction is given by 

Vf is separating C{A) and C{B). We have now to turn toward the indistinguishability conditions for two 
trajectories. Two trajectories are said indistinguishable if they induce the same outputs from the sensor 
network. We have then the following property [2 |. 

Proposition 2. Assume that the sensor network is dense, then two target trajectories (say xt and yt) are 
indistinguishable iff the following conditions hold true: 

r yt = At ±t (At > 0) Vi , 

[ (xt-yt,xt) =0 Vt . 

Proof: First, we shall consider the implications of the indistinguishability. Actually, the two trajectories 
are indistinguishable iff the following condition holds: 

(t, - t„ it) < ^ (t, - U, yt) < Vi V(t„ t,) . (9) 

We then choose tj — ti — a x^ (i.e. ti and tj both belongs to the line separating A and B) and consider 
the following decomposition of the yt vector: 

yt = At Xt + nt , 

so that we have: 

{tj-U,yt) ^a^It\\iU\^ <0. (10) 
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Now, it is always possible to choose a scalar a of the same sign than /i^ So, we conclude that the scalar 
is necessarily equal to zero. Thus , if the trajectories Xt and yt are indistinguishable we have necessarily: 

yt = At xt , . 

Furthermore, the scalar Af is necessarily positive (see eq. |9]l. Then, the lemma [T| inequalities yield: 

{tj-ti,xt) < (xt-yi,Xf) < (tj-tj,Xf). (11) 

Choosing once again tj — = a x^, we deduce from eq. [TT]the second part of prop. |2] i.e. (x^ — y^, Xj) = 
. Considering now the distance between the two indistinguishable trajectories, we have [yt = At Xt) : 

^ ||xt - ytll' = 2 (xt - yt, it - yt) - , (12) 
so that we have Ijxf — yt|l = est. 

Reciprocally, assume that the two conditions yt — \t xt and (xt — yt,xt) = hold true \/t, are the two 
trajectories then indistinguishable? It is sufficient to remark that: 

(yt,yt) = (xf + (yt - xt),yf) = (xt,yt) = A* (xt,xt) , 

(13) 

(ti,yt) = At (ti,xt) . 

Since the scalar At is positive this ends the proof. 

□□□ 

Let us now consider the practical applications of the above general results. 

Rectilinear and uniform motion 

Admitting now that the target motions are rectilinear and uniform (i.e. Xf = xg +^ x). Then prop. |2] yields 
y = A X (A > 0) and: 

(yt - Xt, ±) = (yo - xo, i) + i (1 - A) \\±f = . (14) 

Then, from eq.[T4]we deduce that A = 1 and yo = xo + a x^. So that, the target velocity is fully observable 
while the position is uniquely determined modulo a a x^ translation. 

leg-by-leg trajectory 

Consider now a leg-by-leg trajectory modeling. For a 2-leg one, we have for two indistinguishable 
trajectories: 

" Xt = Xo + ti vi + (t - h) wl , 

where is the velocity of the x(t) trajectory on the i-th leg and ti is the epoch of maneuver. Furthermore, 
we can assume that ti < t^. Considering the implications of prop. |2]both for t < ti and for t > t^, we 
know that if the trajectories are indistinguishable we must have: 

vi = vi and: = . (16) 

So, our objective is now to prove that we have also ti = t^. Considering prop. |2] we thus have the following 
system of equations : 

(yo - Xo + (t - h) (vi - v2) , vi) = for : h < t < t[ (a), 
(yo-xo + (t-ti) (vi-v2),v2) =0 for: h < t < t[ (b), (17) 
(yo - xo + (t'l - h) (vi - v2) , v2) = for : t[ < t (c) . 

Now, on the 1-st leg we have also (yo — xo, vj) =0 (see prop. |2] for t — 0), so that eqs [TTk.b yield: 

(K-v^),vi) = ((vi-v2),vi)=0. (18) 



This means that v]. and are both orthogonal to the same vector (vl. — v^), so they are colHnear, and 
we straightforwardly deduce from eq. [18] that — v^. Finally, it has thus been proved that ti — and 
this reasoning can be extended to any leg number The observability requirements having been considered, 
we turn now toward the development of the algorithmic approaches. Let us first introduce the following 
functional. 

III. The stairwise functional 

Our fist aim is to estimate the target velocity, within a batch processing framework. We assume that 
binary ({ — , +} sensors are uniformly distributed on the field of interest (see fig. |2]l. 
Each sensor will be coupled with a counter, that will be increased by a unity each time-period the sensor 
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Fig. 2. A scenario of target evolution and sensor networic information 

gives us a {+}, an will keep its value each time the sensor gives us a { — }. Then, at the end of the trajectory, 
each sensor has a entire value representing the number of periods the target was approaching. Within a given 
batch, the outputs of the sensor counters can be represented by a stairwise functional (see fig. [Sj. 
Then, once this stair is built, we can define what we call the velocity plane. This plane is the tangent plane 
of the stairwise functional, which means that its direction gives the direction of the stair, while its angle 9 
gives the slope. The direction of the plane gives us the target heading, while the target speed v is given by: 

Thus, estimating the velocity is equivalent to estimating the velocity plane parameters. Mathematical 
justifications are then presented. The target moves with a constant velocity v. Considering the results of 
section |V] its starting position is given by the following equation: 

x(0) = xo + A , A e M , 

so, that : (20) 
x(t) = xo + A v-"- + t V 

This means that at each time period t e M_|_, the possible positions x(t) define a (moving) straight line, 
whose dkection is v-*-. Let us consider now the scalar product (x(/:), v), then we have: 

|^(x(i),v)=«^ (21) 
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This is clearly constant, which means that the surface is a plane. The conclusion follows: the stairwise plane 
is an exhaustive information for the velocity vector We provide in the next section two solutions to estimate 
the velocity plane from the observed data, and give some asymptotic results about the estimation. 

IV. Statistical Methods to Estimate the Velocity Plane 

We showed that estimating the velocity plane allows us to estimate the velocity vector Wile there exists 
several methods to do that, we shall focus on two of them. 

A. The Support Vector Machine ( SVM) approach dSjJ 

As seen previously, the problem we have to face is to optimally separate the two classes of sensors (i.e. 
the + and — ). So, we can use the general framework of SVM, widely used in the classification context. 
The set of labeled patterns {(yi, xi, • • • , y/, x/} (y^ e { — 1, 1} and x^ sensor positions) is said to be linearly 
separable if there exists a vector w and a scalar b such that the following inequalities hold true: 

r (w,x,)+6>l if:y,; = l, 

< (22) 

[ (w,x,>+6<-l if:?/, = -1. 

Let 'H(w, &) = {x|(w,x) + b = 0} (w: normal vector) be this optimal separation plane, and define the 
margin (marg) as the distance of the closest point x^ to H, then it is easily seen that marg — j^^ - Thus, 
maximizing the margin lead to consider the following problem: 

min t(w) = ||w||^ , 

w,b (23) 
s.t. : y, ((w,x,) +6) > 1 Vi = I,- - • y., ±1 . 

Denoting A the vector of Lagrange multipliers, dualization of eq. |23] leads to consider again a quadratic 
problem, but with more explicit constraints [3 1, i.e. : 

max W(A) = -- D K + K^ 1 , 
A 2 (24) 

s.t. : A > , K^Y = , 
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where 1 is a vector made of 1 and Y'^ = (j/i, • • • ,yi) is the /-dimensional vector of labels, and D is the 
Gram matrix: 



(25) 



The dualized problem can be efficiently solved by classical quadratic programming methods. The less-perfect 
case consider the case when data cannot be separated without errors and lead to replace the constraints of 
eq. |23]by the following ones: 



((W,X,) +6) >l-^^A^>0, i = l,--- ,1 



(26) 



Consider now a multiperiod extension of the previous analysis. Let us restrict first to a two-period analysis, 
we shall consider two separating hyperplanes (say Hi,H2) defined by: 



(w, xl) + 6i ^ ±ci according to: yj = ±1 , 
(w, xf) + 62 ^ ±C2 according to: yf — ±1 . 



(27) 



It is also assumed that these two separating planes are associated with time periods T and T + AT, AT 
known. It is easily seen that the margin for the separating plane Hi is while for the plane ^2 it is 
11^. Thus, the problem we have to solve reads: 



mm 

W,Ci,C2,hl,&2 



max 

1:2 



(28) 



s.t.: 



yf ((w, xj) + 61) > ci , yf ((w, xf) + 62) > C2 VZ. 



At a first glance, this problem appears as very complicated. But, without restricting generality, we can 

II w 1 1 

- — ^ — . Making the changes — w — > w 

„2 o ci 



assume that c\ < 02- This means that max 

1,2 

and — ^ bi then leads to consider the classical problem: 



v' .b[,b'2 



S.t. 



yi 



, Xi, 



((w',: 



> 1 VL 



(29) 



Let w* be the (unique) solution of eq. |29] then a straightforward calculation yields the distance d{Hl,H2) 
between the two separating planes, i.e.: 



\bi 



Finally, we deduce that the estimated velocity vector v is given by: 

v = aw* and: ^ = , H^) . 



(30) 



The previous analysis can be easily extended to an arbitrary number of periods, as long as the target trajectory 
remains rectilinear Another definite advantage is that it can be easily extended to multitarget tracking. 

1) 3D-SVM: We can also mix the SVM ideas with that of section Hill Indeed, instead of focusing on a 
2-D dataset, we can consider a 3-dimensional dataset (sensor coordinates and values of the sensor counters). 
The second 3-D dataset is the same, but the value of the counter is increased with unity. So, the separation 
plane is 2-D, and will be as closed to the velocity plane as the sensor number can allow. See fig. |4]for a 
more explicit understanding. The results of the SVM estimation of the velocity plane are discussed in the 
Simulation Results section. 
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Fig. 4. The theoretical stairway of the trajectory. 



B. Projection Pursuit Regression 

The projection pursuit methods have first been introduced by Friedman and Tuckey |4|. Then, they have 
been developed for regression with the projection pursuit regression (PPR) by Friedman and Stuetzle ||5|- 
PPR is mainly a non-parametric method to estimate a regression, with however a certain particularity. Indeed, 
instead of estimating a function / such as Yi — f{Xi) + Ei, where Xi and Yi are known, and Si assuming 
to follow a certain law, PPR estimates g such as Yi — g{Xi 6) + e^. The first step of the algorithm is to 
estimate the direction 9, and then g. In our specific case, 9 will represent the direction of the target, and g 
will give us the value of the velocity. 

1) Modeling : Let Yi be the value of the z-th sensor counter. Xi are the sensor coordinates. If n{Xi 9) 
is the value of the counter i at the end of the track, and p the probability to have the right {+, — } decision, 
we then have {B: binomial): 

C{Y,\X,9)=B{niX,9),p) (31) 

Assuming in a first time that p = 1, the two parameters we would like to estimate are the 9 parameter and 
the n(.) function. 

2) The PPR method in the network context: We have some additional constraints on n{.). First of all, it 
only takes integer values. Then, it is an increasing function (because p = 1). The optimization problem we 
have to solve is the following: 

9 = argmm ^ {n{X, 9) - Y,f , (32) 

where h is calculated in a quite special way. First, we define a non parametric estimation of a function f, 
via: 

\^Kh(Xi 9-u) 

Then, we sort {X 9)i into a vector from the smallest to the biggest. After which we define h{.) via: 

= if /(X >/(X^?(,_i)), ^34^ 

n{X 9(^i-j) = f{X otherwise . 
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Sometimes, due to the integer value of the estimated n(.) function, we have to deal with many possible 
values of 9. Then, in this case, we choose the mean value of 9. Due to the specific behavior of our target 
and our modeling, we know in addition that the general form of n (say n) is given by: 

= ^i\xe)^+(i-i)v,(x9)^+iv]{u) ■ (35) 
The next step is then to estimate v. Such an estimation is given by the following optimization program: 

V = axgmmy^ {n{X9i)-Yif (36) 

3) Convergence: We will study if the estimation is good with an infinite number N of sensors. Assuming 
we have an infinite number of sensors in a closed space, this means that each point of the space gives us 
an information {+,—}. We then will have the exact parameters of the stairwise functional. To that aim, we 
will show in the following paragraph that the probability of having a sensor arbitrary close to the Umits of 
each stair steps is 1. We assume that the sensor positions are randomly distributed, following an uniform 
law. Then, y being fixed: 

C{X\Z) = ^B.„,;S3„,] (37) 

If the velocity vector v is denoted with [a; 6], then: 

Binf = -fy-^, Bsup = (38) 

where (cinf , Cgup) only depends on v and xq, which means that they are deterministic, and independent from 
X. It is quite obvious that Bi^f represents the smaller x-limit of a step, when i?sup represents its higher 
x-hmit. Then, considering the velocity plane, -Binf and -Bgup both belong to the plane. Denote u = inf 
then: 

Ve > P{\u - Binf I < £) = P{u - Binf < e) , (39) 

= P{u<e + Binf) ■ 

where we note A= \u — Binf | < £• We know that: 

r 

P{mfXi <t) = l 1 - if t e [Binf; B,np] , (40) 




Then, we have the following probability calculations: 

P{A) = P{\U - Binf I < s) , 

= P{u<e + Binf) , 

1 ^ -^sup ~ (g + -Binf) ^^Af-, „ XX 

= 1 - ( 5 ZTr > l[B,„f ;-B.,up] ((£ + ^inf )) , 

-'^sup -Oinf 

if e < 

(1 - B.J-B^J '' if ^ e]0; B,np - Binf] 

if e > Bgup — Binf • 

Given the above equation, 1 — -g — 3^-7 is smaller than one, which means that (1 — -g — ^ )^ converges 
to as A'' increases to infinity. Thus, we have finally: 

Ve>0 lim P(|u - Binf | < e) = 1- (41) 

ending the proof. 
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V. Non-linear trajectory estimation 

A. Target Motion Model 

The target is assumed to evolve with a Markov motion, given by: 

Xk|xk-i - AA(FfcXk-i,Qk) (42) 

for fc = 1, 2... where 7V^(/i, ct^) is a gaussian distribution with mean /x and variance a^. The starting position 
is assumed to be unknown. 

B. Sensor Measurement Model and Analysis 

At each time period, each sensor gives us a {+, — } information, meaning that the target is getting closer 
or moving away. Given all the sensors reports at the time-period t, we can easily define a space where the 
target is assumed to be at this time-period. This is the fundamental uncertainty we have at t, and the area 
of this domain is, of course, directly related to the network parameters (sensor number, network geometry, 
etc.). 

C. Velocity Estimation 

We can estimate the direction of the target based on the simple information given by the sensors. Obviously, 
that estimator will only be precise if the number of sensors is significantly great. To perform that estimation, 
we can use several methods, such as the Projection Pursuit Regression Method, or the Support Vector 
Machine Method. The SVM method chosen for our algorithm as a most common method, and is presented 
in the next paragraphs. 

1 ) The effect of target acceleration: To illustrate the effect of velocity change for estimating the target 
position, let us consider a very simple example. Assume that the target motion is uniformly accelerated, i.e. 

Xt = xq + i xo + xo . (43) 

We have now to deal with the following question: Is the target trajectory fully observable? To that aim, 
we first recall the following result. Considering a dense binary network, two target trajectories are said 
indistinguishable iff they provide the same (binary) information which is equivalent to the following 
conditions: 

{ it=yt,(yi-xt,yt)=0 Vt . (44) 

Expliciting the second condition of eq |44l with the target motion model |43] we obtain that the following 
condition holds (Vt): 

(yo -xo,yo> +t{yo -xo,yo) + ^t'^{yo -xo,yo) , 

+t{yo - xo,yo> + t'^{yo - xo,yo) + ^t^ivo - xo,yo> = . 

Thus, (yt — X(,xt) is a zero polynomial, which means that all its coefficients are zero. For the 
coefficients we obtain the condition (yo — XQ.yo) = 0. Similarly with the {yt — xt,Xf) — 0, we obtain 
(yo — xo,xo) = 0. Subtracting these two equalities yield ||yo — xo|| = , or xo = yo- 

Quite similarly, we obtain the equality xo = yo and the last equality: 

(yo -xo,yo +tyo> = . (46) 

Assuming that the couple {yo,yo} spans the sensor space then we deduce that xq = yo. So, it has been 
shown that it was the target acceleration which render the problem fully observable. This reasoning can be 
extended to a wide variety of target modeling. 
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D. Tracking algorithm 

The main issue with the SVM estimation is that it only provides us the general direction of the target 
within a deterministic framework. Moreover, it is highly desirable to develop a reliable algorithm for target 
tracking (velocity and position). To solve this problem, we build a two-step algorithm. In the first step, we 
perform a correction through the estimated unitary velocity vector at each time-period t, called At. Then, in 
a second time, we perform a correction through the orthogonal-estimated (unitary) velocity vector, also at 
each time-period, called Ot- These two corrections give us a better estimation of both the velocity and the 
position of the target. We refer to fig. |5]for the presentation of the rationale of the two correction factors. 

1) The X factor: To build that correction factor, we started with a very simple assumption. At each period 
t, the sensors provide binary motion information. Thanks to the first part of this article, we know that the 
target is in the (special) set lying between the two same-sign-sensors set. Then, starting from the previous 
estimated position of the target, we move the estimated target through the estimated velocity vector direction 
until it stands in that special set. We now define this operator in a mathematical way: 

Let Vt the estimated normalized velocity vector at time t. 

Moreover, let {t| (respectively {t,-^'}j) the coordinates of the sensors (s^) giving a {— } (respectively 
a {+}) at time t. 

We sort vs\ = (vt,t| •*) (respectively wSj-^-* = (v^jt^^')). Then, following a very simple geometrical 
reasoning, we note that (vt;Xt) should be between vsmlx and vs^^}^-^. To ensure that property, we define 
the following correction factor: 

, iisj^jjj"' — (vt,xt-i) 

~ {vt,vt-l> ' 

with the following definition of vsmoy^ ■ ^^^^ 

To calculate this factor, we consider the projection equality: 

(vt, (xt_i + At vt_i)) = vs\+^y^ (48) 

which means that the projection of the corrected value is equal to the mean value of the projection. 
Geometrically, this means that the position of the target is estimated to be in the center of the special 
set defined by the sensors. The value of the correction factor At (see eq. WT\ is then straightforwardly 
deduced from eq. |48] Similarly, the target position is updated via: 

xr"" - + At vt-i . (49) 

Here the correction factor At has been calculated via the average value of the projection. This is an arbitrary 
choice and we can consider the lower or the upper bound of the projection with no significant difference on 
the results of the algorithm. 

Obviously, if the estimation of the position is not very good, the estimated velocity value (clearly based on 
At) will be quite different from the real value of the velocity. The next correction factor is based on the 
assumption that the target velocity changes are upper and lower bounded. 

2) The correction factor: We assume that the velocity of the target has bounded acceleration. Then, if 
the velocity estimated at a certain time t is too different from the velocity estimated at time t—1, this means 
that the estimated position of the target is far from the right one. Then, in that precise case, we consider an 
orthogonal correction, through vt^. 

For that deterministic algorithm we decided to perform a very simple modeling of the velocity. Indeed, we 
take as a right value for the velocity the simple mean of the k previous values of the estimated velocity 
(n^t,k)- We calculate in addition the variance (at _fc), and the factor 9t can be non-zero iff the estimated value 
of the velocity at time t is not in the interval given by [rrit k ^ '^t,k', ^t,k + o't,fc]- We then look for 9t such 
that: 

(x^'' + Ot v,^ - (xt_i + 9t vt^--^; vt-i) - mt,fc . (50) 

The previous equation needs some explanation. Given that Xf is the estimated target position at time t, 
we would like to correct the value to be closer to the right position. The only way we can deal with it. 
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Fig. 5. Correction scenario. 



is to correct the estimated value of the velocity, x™'''' — xt_i is the previous calculated correction. If the 
difference between that estimation and the value mt.k is too important, we try to reduce that difference with 
a translation of the positions at time periods t and i — 1. As we want the positions to stay in the special set 
defined by the sensors, the direction of that translation is given by for the position at time t, and v^j^ 
for the position at time t — 1. 
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Performing straightforward calculation, leads to consider the following correction factor: 

= . (51) 

Obviously, as we could expect when presenting the method, if the target motion is rectihnear and uniform 
, no correction factor can be calculated. Then, the final estimated position is given by: 

^f^n ^ ^corr ^ ^± (52) 

3) The final correction step: Noticeably the most important step of the algorithm, i.e. the 6 correction 
factor, is based on the estimation of the velocity change. Indeed, the best the estimation of the velocity is, 
the best we can estimate the position. Then, our aim is to perform a better analysis of the target motion. 
Considering that from time to time, the estimation of the position increases in quality, a promising way 
should be to perform a feedback of the newest corrector to the oldest position estimation. We denote zt the 
updated estimated position of the target at time t. Then, according to the previous paragraph, the estimated 
position is updated via: 

t 

Vj < t : z, = xf " + e,wi . (53) 

With this new estimator we will be able to perform a better analysis of the target motion (position and 
velocity). 

4) The final algorithm: With the definition of the correction factors, the theoretical part of the algorithm 
is finished. Then, it is presented as follows, at time period t: 



1) Get the binary information of each sensor, and then the target position set. 

2) Estimate the velocity direction at time t via a SVM method 

3) Perform the A calculation, and add that correction to the estimated velocity at time t — 1. The 
time-<-position is then updated. 

4) Check if the estimated velocity at time < — 1 is too different from the modeled value, and in this case, 
calculate Ot- 

5) Update the position at time t, and in this case, the velocity at time t — \ with the correction 9. 
Step 2 and 3 can be inverted with no damage in the process. This is the main part of the algorithm. However, 
there is no mention in that enumeration of the initialization. There are two main state vectors that have to 
be initialized. The position and the velocity. The position is assumed to be unknown, but thanks to the 
sensors, we can have a space where the target is assumed to be at first. We use here a uniform law for the 
initialization, given that we have no further information about where the target can start. 

The initialization of the velocity is not far from that solution. Indeed, with the binary information, we can 
provide a convenient estimate of the velocity direction. Even if we don't have a precise idea of the speed 
value, we can then start the algorithm. 

VI. Simulation Results 

A. Constant Velocity movement 

We shall now investigate the previous developments via simulations. The first figure (fig|6) will show the 
stair built by the previously explained method {N: fixed). The position of the sensors are considered random, 
following a uniform law on the surveillance set. To evaluate the performance of our methods, we decided 
to calculate the mean square error of the two estimated parameters, which are the velocity value and the 
velocity direction. Fig. [T] shows the two MSEs values for both direction and velocity values, assuming the 
sensor number is growing from 10 to 100, and the velocity vector is the [1, 2] vector (m/s). 
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Fig. 6. 2D sensor position and velocity's direction projection. N=100 sensors. Velocity' Vector is [1,2] 
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Fig. 7. Mean Square Error of the Velocity Estimators. Green for the SVM, Blue for the PPR. 



Providing 2000 simulations, the MSEs seems to be unstable. However, the two parameter estimation methods 
leads to a very different conclusion. In the case of the direction estimation, the PPR method works highly 
better than the SVM method, and seems quite stable as the sensor number N grows. On the other side, 
the SVM method is more erratic. One possible explanation is that the PPR method has been first developed 
for the particular case of direction estimation, while the SVM method is more focused on the margins 
maximization, which means in our case a simultaneous estimation of both parameters. The conclusions we 
can make on the velocity value estimation are rather opposite. The MSE becomes reasonable only for the 



14 



SVM method, and for a number of sensors up to 60. Indeed, we have a 0.05 m/s error on a velocity value 
estimation for a theoretical value of a/5- As erratic as the SVM's MSE was in the direction estimation, it 
was however less erratic than the result we have for the PPR value. 

One answer to the MSE erratic value for the PPR could be to find a best way to estimate the velocity value. 
Indeed, in our case, we choose for estimating functional a sum of indicators functions. However, it is not 
clear that this optimization gives a single minimum solution. There could be a finest functional that could 
lead to a most robust optimization solution, and this would be the subject of future works. 

B. Random walk 

We will present in that section the results of the tracking algorithm. We consider here that the target starts 
from the [100, 100] position and that its initial velocity vector is the [1, 1] vector The number of sensors is 
equal to 70, in a quite wide space (300mx300m). The variance of the target motion is not very important, 
and the tracking duration is T = 30 seconds. 

One simulation is presented in figure |8] In red is represented the real target trajectory, quite diffusive, and 
in green the estimated successive positions. The initialization is not very bad because the number of sensors 
is quite important, which means the uniform set is not too large. After the first step, the estimation seems 
to hang the real trajectory, and follows the target well (less than 10 meter error). However, when the target 
turns right, we loose some precision, mainly because the correction factors seemed to be "lost". The reason 
for that behavior is that the SVM method provides us a bad estimation of the velocity vector Then, the 
algorithm provides a correction in a bad direction, which moves away from the real trajectory. During a 
few seconds, the estimation works quite bad, before hanging again the target direction, and then performing 
a quite good estimation of the velocity. Unfortunately, there is no evidence in that example that increasing 
indefinitely the tracking duration results in an estimated position closer and closer to the real target position. 
This is precisely the aim of the two next figures in |9] The first one shows the mean square error of the 



o o o o o o o 




o o o o o o o 



Fig. 8. Trajectory Estimation of a target. In Red, the real trajectory, in green the estimated one. 

estimated position of the target through the trajectory. The total time is T = 30 seconds, and we can see an 
amazing and remarkable decrease of that MSE in the first seconds. It seems however that there is a limit 
to that decrease. Indeed, the MSE will not converge to a zero value, even if we could perform a long-time 
tracking. Clearly, the limitation is due to the binary information at first, and certainly to the number of 
sensors in a second time. Some further work could certainly exhibits a strong link between the number of 
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sensors and the MSE of the position. 

In the same way, the velocity estimation has some acceptable MSE through the tracking process. Despite 
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Fig. 9. MSE of the estimator of the target position. MSE of the estimator of the target velocity. 

the clearly strong decrease at the beginning, the curve then stands to an acceptable but non zero value. The 
effect is more obvious than in the position case, surely because of the velocity modelling we make in the 
algorithm, which forces the velocity estimation to very bad evolution. A clue could be to perform a most 
sophisticated modelling of the velocity, but given the binary information, this won't be easy. This is another 
work in progress for the evolution of our algorithm. 

VII. Conclusion 

A. CVM 

In this paper, we chose to focus on the use of the { — , +} at the level of information processing for a 
sensor network. Though this information is rather poor, it has been shown that it can provide very interesting 
results about the target velocity estimation. The theoretical aspects of our methods have been thoroughly 
investigated, and it has been shown that the PPR method leads to the right velocity plane if the number of 
sensors increase to infinity. The feasibility of the new concept ("velocity plane") for estimating the target 
trajectory parameters has been put in evidence. The proposed methods seem to be sufficiently general and 
versatile to explore numerous extensions like: target tracking and dealing with multiple targets within the 
same binary context. 

B. Random walk 

A new method for tracking both position and velocity of a moving target via binary data has been 
developed. Though the instantaneous data are poorly informative, our algorithm takes benefit of the network 
extent and density via specific spatio-temporal analysis. This is remarkable since the assumptions we made 
about target motion are not restrictive. Noticeably also, our algorithm is quite fast and reliable. Furthermore, 
it is clear that performance can be greatly improved if we can consider that the acquisition frequency is 
(far) greater than the maneuver frequency. In particular, we can mix the present method with the one we 
developed in |[l]. 
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However, some important questions remain. The first one concerns the velocity modeling. We focused 
on this paper on the adaptability of the different correction factors, but we didn't pay much attention to 
that modeling, which can definitely improve the estimation quality. Moreover, our tracking algorithm is 
basically deterministic even if the target motion modeling is basically probabilistic. Thus, it should be worth 
to calculate the first correction factor (A) via a likelihood, such that x'^"'"'" does not always stays in the 
mean of the special set. Moreover, that likelihood should be related to all the sources of sensor uncertainty. 
In addition, the present algorithm gives a slow response to sudden target maneuver. A remedy should be to 
incorporate a stochastic modeling of such event in our algorithm. 

The second correction factor (6) may also be improved via a stochastic approach. Instead of considering 
a correction only related to the estimated velocity estimated, we could immerse this correction within a 
stochastic framework involving both and 6. These observations are part of our next work on that very 
constrained but also quite exciting tracking framework. The last important point is multiple target tracking. 
Even if our work in this area is quite preliminary, it is our strong belief that our spatio-temporal separation 
based algorithm should be the natural way to overcome the association problems. 
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